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Abstract We study the properties of galaxies at redshift z — 2 in a A cold dark matter 
(ACDM) universe, using two different types of hydrodynamic simulation meth- 
ods - Eulerian TVD and smoothed particle hydrodynamics (SPH) - and a spec- 
trophotometric analysis in the U n ,G,R filter set. The simulated galaxies at z = 
2 satisfy the color-selection criteria proposed by Adelberger et al. (2004) and 
Steidel et al, (2004) when we assume Calzetti extinction with E(B— V) = 0.15. 
We find that the number density of simulated galaxies brighter than R < 25.5 at 
z = 2 is about 1 x 10" 2 h 3 Mpc" 3 for E(B - V) = 0.15, which is roughly 
twice that of the number density found by Erb et al. (2004) for the ultravio- 
let (UV) bright sample. This suggests that roughly half of the massive galaxies 
with M* > 10 10 h~ 1 M e at z = 2 are UV bright population, and the other 
half is bright in the infra-red (IR) wavelengths. The most massive galaxies at 
z — 2 have stellar masses > 1O 11_12 M0. They typically have been contin- 
uously forming stars with a rate exceeding 30 Mq yr _1 over a few Gyrs from 
z = 10 to z = 2, together with significant contribution by starbursts reaching up 
to 1000 Mq yr _1 which lie on top of the continuous component. TVD simula- 
tions indicate a more sporadic star formation history than the SPH simulations. 
Our results do not imply that hierarchical galaxy formation fails to account for 
the observed massive galaxies at z > 1. The global star formation rate density 
in our simulations peaks at z > 5, a much higher redshift than predicted by the 
semianalytic models. This star formation history suggests early build-up of the 
stellar mass density, and predicts that 70 (50, 30)% of the total stellar mass at 
z — had already been formed by z = 1 (2, 3). Upcoming observations by 
Spitzer and Swift might help to better constrain the star formation history at high 
redshift. 



1. Introduction 



A number of recent observational studies have revealed a new population of 
red, massive galaxies at redshift z ~ 2 (e.g. Chen et al. 2003; Daddi et al. 
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Figure 1. Stellar mass vs. G — R color of simulated galaxies at z = 2. Gray points show the 
total galaxy population, and the black crosses indicate the galaxies brighter than R = 25.5. The 
vertical long-dashed lines roughly indicate the UV color selection range of Steidel et al. (2004). 



2004; Franx et al. 2003; Glazebrook et al. 2004; McCarthy et al. 2004; van 
Dokkum et al. 2004), utilizing near-IR wavelengths which are less affected by 
dust extinction. At the same time, some studies focused on the assembly of 
stellar mass density at high redshift by comparing the observational data and 
semi-analytic models of galaxy formation (e.g. Dickinson et al. 2003; Fontana 
et al. 2003; Poli et al. 2003). These works have suggested that the hierarchical 
structure formation theory may have difficulty in accounting for sufficient early 
star formation. The concern grew with the mounting evidence for high redshift 
galaxy formation including the discovery of Extremely Red Objects (EROs) 
at z > 1, sub-millimeter galaxies at z > 2, Lyman break galaxies (LBGs) at 
z > 3, and Lyman-a emitters at z > 4. We now face the important question as 
to whether this evidence for high-redshift galaxy formation is consistent with 
the concordance ACDM model. See Nagamine et al. (2004a,b) for the details 
of the simulations and the present work. 

2. Results 

Figure 1 shows the stellar mass vs. G-R color of the simulated galaxies at 
z = 2. The most massive galaxies have stellar masses M* > 10 10 /i _1 M Q , and 
UV colors — 1.2 < G— i? < 0.8, consistent with the UV color selection criteria 
of Steidel et al. (2004). The differences in the distribution of the points can be 
understood in terms of the different box sizes and the randomness of the initial 
condition of the simulations. On the right-hand-side of the panels, there are 
a couple of red, passive systems that have stellar masses M* > 10 10 /i _1 M Q . 



Massive galaxies at z = 2 
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Figure 2. Star formation history of the most massive galaxies (panels (a) & (c)) and the 
reddest galaxies (panels (b) & (d)) in the simulation. The left two panels are for the SPH G6 
run, and the right ones are for the TVD run. 



The near-IR properties of these passive systems will be reported in future pa- 
pers. Figure 2 shows the star formation histories of the most massive and 
most reddest galaxies in the simulations. For the most massive systems, the 
star formation rate has a continuous component of > 30 M Q yr _1 over a few 
Gyrs, and a starburst component exists on top of the continuous component 
that reaches up to 1000 M Q yr _1 . Such extreme star formation histories al- 
low these galaxies to build up stellar masses larger than a few xlO 11 M Q by 
z = 2. The TVD simulation indicates somewhat more sporadic star forma- 
tion history, which is perhaps due to the differences in the details of the star 
formation recipe and the resolution. 

Figure 3 shows the star formation history of the entire simulation box. In 
panel (a), all the models including the analytic model of Hernquist & Springel 
(2003) show that the SFR density peaks at z > 5. These SFR histories lead to 
early build-up of the stellar mass density compared to both the current obser- 
vational estimates and the results of the semianalytic models (see Nagamine 
et al. 2004b for a direct comparison), and we predict that 70 (50, 30)% of the 
total stellar mass at z = had already been formed by z = 1 (2, 3) based on 
these theoretical models. 

In summary, we have shown that the simulations based on the hierarchi- 
cal CDM model can in fact account for the masses and the comoving number 
densities of the massive galaxies at z = 2 that are found by the recent observa- 
tions. Our simulations indicate that the properties (i.e. stellar mass, color, SF 
history, clustering) of the UV bright LBGs at z > 2 can be understood if they 
are identified with the galaxies that reside in the massive dark matter halos at 
the time. 
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Figure 3. Panel (a): Star formation rate density as a function of redshift for simulations and 
the analytic model by Hernquist & Springel (2003). The sources of the extinction corrected data 
points are described in Nagamine et al. (2004a). Panel (b): Stellar mass density as a function of 
redshift. 
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